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ABSTRACT 

The material that is ejected in a common-envelope (CE) phase in a close binary system 
provides an ideal environment for dust formation. By constructing a simple toy model to describe 
the evolution of the density and the temperature of CE ejecta and using the AGBDUST code to 
model dust formation, we show that dust can form efficiently in this environment. The actual dust 
masses produced in the CE ejecta depend strongly on their temperature and density evolution. 
We estimate the total dust masses produced by CE evolution by means of a population synthesis 
code and show that, compared to dust production in AGB stars, the dust produced in CE ejecta 
may be quite significant and could even dominate under certain circumstances. 



Subject headings: binaries: close — stars: evolution 

1. Introduction 

Dust is one of the important ingredients of 
the interstellar medium (ISM). It plays a central 
role in the astrophysics of the ISM, from the 
thermodynamics and chemistry of the gas to the 
dynamics of star formation. It affects the thermal 
and chemical balance of the ISM by reprocessing 
the radiative output from stars, providing photo- 
electrons which heat gas, and depleting the gas of 
refractory elements, whi ch are important co oling 



agents in the ISM (e.g.. iNozawa et al.l 120011) . In 

addition, dust changes the spectra of galaxies: 
radiation at short wavelengths is attenuated, an d 



energy is radiated in the infrared. iBernstein et al 



(|2002l ) estimated that 30% or more of the energy 
emitted as starlight in the Universe is re-radiated 
by dust in the infrared. 

From a theoretical point of view, the formation 
and growth of dust grains is still a wide l y unsolved 
probl em (e.g., Gail &: Sedlmavr 19991 : Gall et al. 
20111 ). There are numerous efforts under way 



trying to underst and the processes involved(e.g. 
Gail fc Sedlmavd Il999l : iTodini fc Ferraral 12001 
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circumstellar matter — dust 



Draine l2009l) . Up to now, people have assumed 
that dust mainly originates from the stellar 
winds of asymptotic giant branch (AGB) stars 
and supernova (SN) ejecta. Gail and his 
collaborators investigated the formation and 
growth o f dust gra ins pr o duced by AGB stars (c.i 



Gail fc Sedlmavn Il999l: 



Gail et al 



200 



Ferrarotti fc Gai l 2p0e 
while iTodini fc Ferraral f20o"l 



Noza wa et al.l (|2003l ) and iBianchi fc Schneider 
(2Q0J7 

among others, have studied the conden- 
sation and survival of dust grains in SN ejecta. 

Common envelopes (CE) form as a result of 
dynamical timescale mass exchange in close bina- 
ries and play an essentia l role in their evolu tion 
(see, e.g., Paczynsk] 1976 : Iben &: Livio 1993 ). In 
most cases, CE evolution involves a giant star 
transferring matter to a normal or a degenerate 
star on a dynamical timescale. The giant envelope 
overfills the Roche lobes of both stars and engulfs 
the giant core and its companion. During the 
CE phase, orbital energy is transferred to the 
CE via dynamical friction between the orbiting 
components and the non-corotating CE. If enough 
orbital energy is deposited in the CE before the 
components merge, the whole envelope can be 
ejected on a dynamical timescale. The CE ejecta 
then rapidly expand. According to th e classical 
theory of nucleation ( Feder et al. 19661) , the sat- 
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uration pressure in most cases falls more rapidly 
than the pressure of the vapour, and the vapour 
becomes supersaturated when saturated vapour 
expands adiabatically. The formation of dust 
grains from the gas phase can occur from vapour 
in a supersaturated state. Therefore, CE ejecta 
provide a potentially important environment for 
the formation of dust grains. To our knowledge, 
there is no investigation of dust formation in CE 
ejecta to date. 

In this paper, we investigate dust grain for- 
mation and growth in CE ejecta. §2 describes 
the model for the CE ejecta, and §3 gives the 
input parameters for the dust modeling. The main 
results are presented in §4 and discussed in detail 
in §5. In §6 we estimate the total dust masses 
produced using population synthesis modeling, 
and §7 summarizes our main conclusions. 

2. The Mass Density and Temperature of 
Common Envelope Ejecta 

In general, a CE system contains a gainer 
(a normal star or compact object) and a donor 
composed of a giant envelope and a giant core. 
The CE ejecta are formed from the ejected giant 
envelope. In general, dust formation and growth 
is determined by the density and temperature 
evolution of the CE ejecta. Unfortunately, our un- 
derstanding of CE evolution is still rather poor de- 
spite num erous efforts to improve our understand- 
ing (e.g.. iRicker fc Taa3 l2008t ICe et al.l I2OIOI: 
Delove fc TaamI 2010 ). To first-order approxima- 
tion, we assume that the whole giant envelope is 
ejected as CE ejecta on a dyn amical timescale. 
We use the EZ code (jPaxto n 2004} derived 



from Eggleton's STARS code (lEggletonlll97lD . to 
simulate the giant's structure and evolution. In 
this work, 'giant envelope' refers to the region in 
the star in which the hydrogen abundance (by 
mass) is larger than 0.5. In the EZ code, the giant 
envelope is divided into ^ 100 zones. The mass 
density and the temperature in every zone are 
noted by po and Tq, respectively. These determine 
the initial conditions for the CE ejecta. 

After the giant envelope has been ejected, 
it begins to expand, and the density and 
temperature start to decrease. 

(i) The Evolution of Mass Density - The mass 



-Md onor 5 

where Mgainer 



is the mass of the 
onor is the mass of the donor. If 



density p is represented hy p — Mcj/(47ri?^y), 
where Mej and R are the mass-ejection rate and 
the radial distance of the ejected matter, re- 
spectively, and V is t he veloc i ty of the ejected 
matter. According to ^Harped ( 1996[ ). the main 
characteristic of the cool winds of evolved K and 
early M giants is that their terminal velocities is 
lower than the surface escape velocity, typically, 

Voo < ^Vcsc = \\f^^^-> where V^sc is the escape 
velocity at the stellar surface, Rs is the stellar 
radius and G is the gravitational constant; Mt = 

-^^gainer 

gainer and Md, 
the velocity of the CE ejecta at infinity is similar 
to the terminal velocity of the cool wind from 
the red giant, energy conservation can be used to 
determine the velocity as a function of distance 
R from the giant from _ y^) = GMt/R. 

Considering the stellar structure of a red giant, 
we assume, for simplicity, that V approximately 

equals \J^^^^-, that is, V is the escape velocity at 

R. It is difficult to determine Mcj for CE ejecta. 
Since we assume that the giant envelope is ejected 
on a dynamical timescale, we assume for each layer 
within the CE ejecta that A/cj = pqAttRqV, where 
Rq is the initial distance of this layer to the stellar 
center. Therefore, the evolution of the density in 
each layer can be represented by 
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(ii) The Evolution of Temperature - The expan- 
sion of the CE ejecta leads to the cooling of the 
gas. At the beginning, most of the hydrogen 
atoms in the giant envelope are fully ionized except 
for hydrogen near the giant's photosphere. As 
the temperature decreases, hydrogen recombines, 
gradually turning ionized hydrogen nuclei into 
hydrogen atoms and releasing the recombination 
energy in the process. Here, to first-order ap- 
proximation, we only consider the recombination 
of hydrogen. The released energy can partly be 
absorbed by the CE eje cta and be used t o drive the 
CE expansion further (jHan et al.lll995l ) or be lost 
from the CE ejecta (e.g. by radiation). However, 
it is difficult to accurately describe the efficiency 
of this process. For simplicity, we introduce a 
parameter 7 to describe the temperature evolution 
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in each layer: 



T = 



R 



To, 



(2) 



where Tq is the initial temperature. The degree of 
ionization of the hydrogen atoms is determined by 
the Saha equation 

n — Uc ^/W2jfmJtBT \kBT J ' 

where n, and no are the density of hydrogen 
atoms, hydrogen ions and electrons, respectively. 
In the case of pure hydrogen, n+ = Ue- The 
constants h, fee and mo are the Planck constant, 
the Boltzmann constant and the mass of the 
electron, respectively. The ionization energy of 
hydrogen is given by e and equals 13.6 eV. The 
temperature T is given by Eq. As the 

temperature decreases, the degree of ionization of 
hydrogen, n^/n, also decreases. In this work, we 
take rij^jn — 1% as the boundary between the 
regions where hydrogen is fully ionized and where 
it is partially ionized or atomic. The distance of 
this boundary from i?o is denoted as i?cr, and the 
temperature of the CE ejecta at i?cr is Once 
most of the ionized hydrogen has recombined, 
the gas in the CE ejecta will be similar to the 
gas at the surfa c e of t he red giant . Fol l owing 
Gail fc Sedlmavil (|l999[ ). we refer to ILucvI (|1971 



19761 ) to calculate the temperature evolution from 
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where tl is defined by 
drL 

Here, kh is the fiux averaged mass extinction coef- 
ficient and is calculated by a simple superposition 
of the extinction of t he different dust species and 
the gas (see details in lGail fc Sedlmavi ( 1999f )). 

In summary, the expansion of the CE ejecta 
is split into two zones. In the inner zone (i? < 
i?cr), the evolution of the CE ejecta's temperature 
is given by Eq. In this zone, the CE 

ejecta temperature is high enough for hydrogen 
to be fully ionized so that dust cannot form. In 
the outer zone (i? > i?cr), the temperature's 
evolution is described by Eq. In this zone, 

the temperature of the CE ejecta has dropped so 
that ultimately dust grains can form. 



3. A Dust Model for the Common Enve- 
lope Ejecta 

As mentioned in the last section, we assume 
that the structure of the CE ejecta in the outer 
zone, where i? > i?cr, is similar to that in the 
stellar wind from an AGB star. With the CE 
ejecta expanding, its temperature and pressure 
decrease, and some tiny seed nuclei can form in the 
cooling gas. However, the nucleation of seed nuclei 
from the gas phase is a difficult problem. The 
AGBDUST code does not consider this problem 
and assumes that the seed nuclei already exist. 
Condensation of dust can occur on the surface 
of the seed nucl ei. Th en, dus t starts to grow. 
[Gail fc Sedlmavi^ ('1999') and iFerrarotti Sz Gaill 
(|2001i . 12002. ,2003 . ,2005) have investigated the 
condensation and growth of dust grains in the 
stellar wind from AGB stars using the AGBDUST 
code. In this work, we use the same code for 
dust formation in CE ejecta. In the AGBDUST 
code, the condensation and growth of dust grains 
is affected by the following input parameters 
(other input parameters that arc not specifically 
m entioned are tak en to have the default values as 
in IFerr arotti fc Gail (200Q)): 

(i) Chemical composition - The chemical com- 
position has a large effect on dust formation. 
In this paper, we investigate dust formation in 
CE ejecta formed in binary systems in which 
giants are on the first giant branch (FGB). Due 
to the first dredge-up, the chemical abundances 
in the envelopes of FGBs star are different from 
the initial abundances, which are taken from 



Anders fc Gre vcssc (1989) for the Sun. The effects 
of the first dredge-up are a reduction of ^^C by 
approximately 30% and no cha nge in the ^^O 



abundance at the stellar surface ( Iben fc Renzinil 



19831 ). For Fe, Mg and Si, which are some of the 



key elements for dust formation, the abundances 
are not changed substantially. Therefore, the 
abundance ratio of carbon to oxygen by number 
(C/0) in th e CE ejecta is approximately 0.4. 
According to ICail fc Sedlmavr (1999), the most 
abundant dust species formed in the circumstellar 
matter (where C/0< 1) are olivine- and pyroxene- 
type silicate grains and iron grains, 
(ii) Temperature ~ The temperature of the CE 
ejecta is determined by Eqs. ([2]) and dJ]) and 
is affected by the parameter 7, a parameter 
that is quite uncertain. Based on the results 
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of model calculations by iFransson fc Chevalier 
iKozasa et al.l I (|l989h adopted an adiabatic 
index 7ad = 1-25 for the early stage of SN 
explosions. For an adiabatically expanding perfect 
gas, pT^-Tad = constant. We assume that the 
temperature evolution of the CE ejecta in the 
inner zone is similar to that in the early stage of 
a SN explosion. This implies 7 — 0.375 in Eq. 
([2]). In this work, in order to check the effects of 
this parameter on the dust-formation process, we 
simulate cases with 7 = 0.2, 0.3 and 0.4. 

(iii) The mass density - The mass density is 
determined by Eq.([T]). The density profile will 
also be affected by the geometry of the CE ejecta. 
[Sandquist et al. 1 (fl998i) showed that the mass loss 
in the orbital plane is about 5 times larger than in 
the polar direction. It is beyond the scope of this 
paper to simulate such more realistic structures. 
Specifically, in this work we assume that the CE 
ejecta expand spherically. 

(iv) The stellar luminosity and mass - The stellar 
luminosity is taken to be the luminosity of the 
FGB star, and the mass is the mass of the 
binary, including the FGB star and the degenerate 
companion. However, their effects are weak 
because the region of dust formation is far away 
from the binary system. 

4. Results 

Using the EZ code, we simulate the evolution 
of 10 stars with initial masses of 1.0, 1.25, 1.5, 
1.7 2.0, 2.5, 3.0, 4.0, 5.0 and 7.0 Mq, respectively. 
Their companions are assumed to be degenerate 
stars of 1.0 Mq. We assume that they experience 
a CE phase immediately when they overflow their 
Roche lobes on the FGB, ejecting their whole 
envelopes. The density and the temperature of 
the CE ejecta depend on the parameter 7, and 
the initial temperature and density profiles Tq and 
Pq. Tq and pa in every layer are given by the EZ 
code and depend on the evolutionary state of the 
giant. In order to examine the effects of Tq and po 
on dust formation, we calculate the dust masses 
produced in the CE ejecta which originate from 
giants at the top of the FGB and at the base of 
the FGB, respectively. Figure [T] shows Tq and po 
at every level of the giant envelope at the top of 
the FGB and at the base of the FGB. 

Figure [2] shows the amount of dust produced (in 



Mq) in the CE ejecta in these models for different 
values of 7. As the figure demonstrates, 7 and the 
evolutionary state of the giant dramatically affect 
the dust masses produced. In particular, varying 
7 from 0.2 to 0.4 introduces an uncertainty for the 
dust masses of up to a factor of ^ 10^. Different Tq 
and Po profiles, i.e. varying the evolutionary state 
from the base of the FGB to the top of the FGB, 
introduce a variation of up to factor of ^ 10^, but 
this depends strongly on the mass of the giant. 

The dust masses mainly depend on the tem- 
perature and the mass density of the CE ejecta. 
In order to show the effects of 7, Tq and po, we 
give the temperature, mass density and relative 
distance of the CE ejecta produced by an FGB 
star of 1.0 Mq at Ra in Figure [3l Obviously, a 
higher 7 results in a smaller _/?„• for the same Tq 
and Po, and the CE ejecta have a higher density 
(see the right p anels or the left panels in Figure 
121) . According tol Cail fc Sedlmavil (|l999[) , a higher 
mass density results in a higher degree of dust 
condensation in the stellar wind. Thus, in the 
simulation with a high 7, the CE ejecta offer a 
very favourable environment for the formation and 
growth of dust grains. Similarly, a lower Tq also 
results in a smaller for the same 7 (see Figure 
[3]), and the CE ejecta have higher density at i?cr. 
Again dust grains can easily form and grow in 
the CE ejecta in this case. In contrast, a low 7 
and a high Tq result in a large Ra-, and the mass 
density of the CE ejecta becomes too low to lead 
to efficient dust production. 

Beyond the CE ejecta enter the zone 

where dust may form. Figure |4] gives the 
quantities of olivine- and pyroxene-type silicate 
grains and iron grains produced by the CE 
ejecta. Due to the small stick ing efficiency for 
quartz (jCail fc Sedlmavd Il999h . its quantity in 
dust grains is negligible. In our simulations, the 
silicates mainly consist of olivine and pyroxene 
whose proportion in the dust grains is between 
70—90 %. The proportion of iron in the dust grains 
is between ~ 10 — 30%. As Figure [4] shows, for 
a small 7 and a high Tq, dust is mainly produced 
by the CE ejecta close to the stellar surface, while 
the whole CE ejecta can efficiently produce dust 
for a high 7, and a low Tq and pq. 

Ferrarotti fc Gail (|2006h calculated the quan- 
tities of dust produced by AGB stars. This 
is plotted as a dot-dashed curve in Figure [2] 
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Fig. 1. — The initial temperature and density profiles of the envelopes of giants (at different evolutionary 
stages as indicated) as a function of relative mass coordinate. 
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This shows that, although dust production is 
comparatively unimportant in the CE models with 
7 = 0.2 at the base of FGB, the dust quantities 
produced by the CE ejecta are significant or may 
even dominate in other models with a higher 7. 
Th ere are two main r e asons for this: 

(i) iFerrarotti fc Gaill (|2006f ) concluded that the 
dust condensation is efficient when AGB stars 
have high mass-loss rates, higher than 2 x 
10~^ Mqjt-^. For an AGB wind with a 
high mass-loss rate of 1 x 10~^Moyr~^, the 
temperature and the mass density of dust-forming 
zones are ~ 1.3 x 10^ K and 5.0 x 10~^^g cm""^, 
respectively. About 32% of Si elements in an 
AGB wind condensate into olivine-, pyroxene- and 
quartz- type dust, and about 4% of Fe elements 
condensate into iron grains. However, in the 
simulation with 7 — 0.4 and at the top of FGB, 
the temperature and the mass density of dust- 
forming zones in the CE ejecta are ~ 1.1 x 10'^ 
K and ^ 3.3 x lO^^^g cm""^, respectively. The 
mass density of the CE ejecta in the dust-forming 
zone is much higher than that in an AGB wind 
and hence is very favorable for dust formation 
and growth. About 90% of the Si elements in 
the CE ejecta condensate into silicate grains, and 
about 50% of the Fe elements in the CE ejecta 
condensate into iron grains. 

(ii) Generally, about 7% (for M = 1.0 M©) 
to 70% (for Mi = 7.0 Mq) of the mass of an 
AGB star are lost at mass-loss rates in excess 
of 2 X 10~^ M0 yr"-'^, and only this portion can 
efficiently produce dust. However, as Figure [1] 
shows, about 40% (for M = l.OM©) to 80% (for 
Mi = 7.0 Mq) of the stellar mass can be ejected 
as a CE at the top of the FGB. Furthermore, 
as the top-left panel in Figure H] shows, dust can 
efficiently form and grow throughout the CE ejecta 
in this case. 

Therefore, dust formation and growth in CE 
ejecta with high 7 is more efficient than in an AGB 
wind, and, under these circumstances, CE ejecta 
may produce more dust than AGB winds. 

In IFerrarotti fc Gaill (|2006l ). due to the third 
dredge-up, the dust species produced by AGB 
stars are olivine-, pyroxene- and quartz- type sil- 
icates, iron, SiC and carbon dust grains. Stars 
with initial masses between ~ 2 — 4 Mq can evolve 
into carbon stars. They produce a large amount 
of carbon dust, which produces the bump in the 



thick dot-dashed curve in Figure |2l In our work, 
because the CE ejecta originate from FGB stars, 
the dust species produced in the CE ejecta are 
olivine-, pyroxene- and quartz- type silicates and 
iron grains. If an AGB star which has undergone 
the third dredge-up experiences a CE phase, the 
C/0 ratio in the CE ejecta may be higher than 1. 
Then, other dust species (such as SiC and carbon 
dust grains) can also form. 

5. Discussion 

As §4 suggests, under certain circumstances, 
CE ejecta can be very efficient producers of 
dust. Hence, it should be possible to observe 
large amounts of dust around post-CE systems. 
Ferrarotti &: Gaill (|2006l ) showed that the distance 
of dust formation from AGB stars is about 
10^^ — 10^'^ cm. Figure [5] shows the amount 
of dust produced by CE ejection as a function 
of distance from the FGB stars. The distance 
of dust formation in the CE ejecta is between 
~ 10^^ — 10^* cm and is relatively far away from 
the FGB star; this may make it difficult to observe 
the dust produced in CE ejecta. 

However, if there is an extreme event after the 
dust grains have formed in the CE ejecta, for 
example, a SN in which a large amount of energy 
is released, dust grains can be illuminated. This 
may produce light echoes, as have been observed 
in several SNe, by radiation scattered by the dust 
near or along the line of sight. 

For example, the progenitors of Type la SNe 
(SNe la) may have experience d CE evolution. 



imagmg. 



Garnavich et al 



From HST WFPC2 ^ ^, 

(|200ll ) proposed that SN 1998bu may have two 
echoes caused by dust at < 10 pc and 120 ± 15 
pc away from the SN: the inner echo is likely to 
be caused by dust from circumstellar material, 
while the outer component is consistent with ISM 
dust. SN 1998bu is a SN la. There are three 
other known SN la echoes, SNe 1991T, 1995E 



and 2006X. IWang et al.l (|2008l) showed that these 
echoes have a wide range of dust distances from 
< 10 pc to ~ 210 pc, which is consistent with 
our results (see Figure [5]). However, there are 
no observational data for light echoes in the 
majority of SNe la. The main reason is that it 
is very difficult to observe these light echoes. In 
the single-degenerate model, there is typically a 
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Fig. 5. — The amounts of dust produced in the CE ejecta as a function of the distance to the FGB star. 
The left panels are for FGB stars with masses less than 2.0 M©, and the right panels are for FGB stars with 
masses larger than 2.0 M©. 
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Fig. 2. — Dust masses as a function of initial 
stellar mass. FG200 6 refers the results in 
Ferrarotti fc Gaill (l2006h . The circles give the dust 
masses in the models calculated. 



long time delay between the CE phase and the 
supernova explosion of betwe e n 10^ and 10^ year s 
(|Han fc Podsiadlowskil l2004t iMeng et all l2009l ). 
In the double-degenerate model, the delay times 
often exceed 10^ years. In contrast, the time it 
takes for material to move from the binary to 
the zone of dust formation is only 10 — 100 yr. 
Dust grains formed in the CE ejecta probably 
have significant outflow speeds. It is difficult 
to determine this outflow speed. If it is ~ 10 
km/s, dust grains will have moved far away from 
the binary by the time of the SN, and it will 
be difficult to observe light echoes. However, 
if a degenerate WDs explode as SNe la within 
10'' — 10* yr after the CE phase, it is possible to 
observe light echoes at distances of 100 — 1000 pc 
(depending on the outflow velocity). The delay 
time in the single-degenerate model depends on 
the mass of the companion of the WD. A delay 
time of ~ 10* yr means that the companions of 
WDs have initial masses larger than ~ 6 Mq. 
Considering the co ndition for dynamical s t ability 
for mass transfer, Han fc Podsiadlowskil ( 2004 ) 
suggested that the main-sequence masses in the 
progenitors of SNe la cannot exceed ~ 4.0 M©. 
However, using observations of the evolution of 
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the SN la rate with redshift, Mannucci et al 



Fig. 3. — Temperature, density and relative 
distance as a function of relative mass coordinate 
of the CE ejecta at Rcr- The left panels are for the 
CE ejecta produced by a star with IM© at the top 
of the FGB. The right panels are for the CE ejecta 
produced by a star with IMq at the base of the 
FGB. Solid and dashed curves represent 7 = 0.4 
and 7 =0.325, respectively. 
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Fig. 4. — Similar to Figure |31 but for the amounts 
of olivine- and pyroxene-type silicate grains and 
iron grains produced as a function of relative mass 
coordinate. 




(|2006[ ) suggested that SNe la have a wide range 
of delay times, from < 10® to > 10^° years. 
From a theoretical point of view, ass uming an 



asphe rical stellar wind in symbiotic stars. iLii et al 



( 2009t ) claimed that the initial masses of WD's 
companions in pr ogenitors of SNe la can exceed 
6.0 M©. Similarly. IWang et al.l (|2009l ) argued that 
the delay times of some SNe la from the helium- 
star channel are shorter than 10® years. Therefore, 
it may be possible in principle to observe light 
echoes in SNe la in which their progenitors are 
composed of WDs and massive companions. 

In one of the main channels for Type lb SNe 
(core collapse supernovae without hydrogen), the 
progenitors are believed to lose their hydrogen 
envelopes in a CE phase ( Podsiadlowski et al.l 
1992i) . The resulting naked helium stars (which 
may appear as Wolf-Rayet stars) will explode after 
the CE ejection within ^ 10^ — 10^ yr (depending 
on the evolutionary state at the beginning of 
mass transfer). Therefore, it may be possible 
to observe light echos from dust grains in SNe 
lb, although their luminosities are much lower 
than those in SNe la. In fact, dust grains have 
been observed in the case of SN2006jc, which 
had a Wolf-Rayet star progenitor, though, in this 
case, it may be more likely that t he dust was 
produce d in the SN ej ecta themselves ('Smith et alj 
2008; P i Carlo et aD 2008; Nozawa ct al. 200^ 
Up to now, to our knowledge, there is no direct 
observational evidence for dust grains in SNe lb 
formed by CE ejection. 

6. An Estimate of the Dust Produced by 
CE Ejecta via Population Synthesis 

CE evolution often occurs in close binary sys- 
tems, and the CE ejecta may provide an im- 
portant contribution to the dust in the ISM. Its 
importance to the overall dust production can 
be estimated using population synthesis model- 
ing. Similar to the main case considered in our 
study of symb iotic stars with WD components 
(iLii et al.ll200(jl. we use the initial mass-function 
of iMiller &: Scald (|l979l ) for the mass of the pri- 
mary c omponents and a flat distribution of mass 
ratios ( Kraicheva et al. 19891 : Goldberg fc Mazeh 
19941) . The distribution of separations is deter- 
mined by log a = bX + 1, where X is a random 



variable uniformly distributed in the range [0,1] 
and the separation a is in i?©. 

During the CE phase, the binary experiences a 
dynamical spiral-in. It is generally assumed that 
the orbital energy that is released by the spiral- 
in process is used to expel the envelope of the 
donor with an efficiency ace- In the theoretical 
calculation, the dynamical spiral-in is affected 
by the combined parameter accAcc, where Aco 
parameterizes the envelope structure of the donor. 
In this work, we take acoAco = 1-0, i.e. assume 
that the CE ejection process is very efficient. 

Using the rapid bi n ary-st ar evolution (BSE) 
code of Hurlev et al.l ( 2002 ). we calculate the 
evolution of 10^ binary systems. About 20% of 
all binary systems undergo CE evolution and eject 
their envelopes. Of these, about 51 % come from 
FGB stars, while the rest originate from stars 
on the AGE. For simplicity, we assume that all 
ejected CEs originate from the stellar envelope 
either at the base of the FGB or the top of 
the FGB. Usually, the envelopes of AGB stars 
have temperature lower than those of FGB stars. 
Therefore, this assumption may underestimate the 
amount of dust produced. According to Figure 
[21 for given input parameters, the dust masses 
mainly depend on the stellar mass at the beginning 
of the CE evolution. Using one-dimensional linear 
interpolation of stellar masses, we estimate the 
quantities of dust produced by CE evolution using 
this population sy nthesis approach. Simi larly, 
using the results of iFerrarotti fc Gaill (l2006l) . we 
also estimate the quantities of dust produced by 
AGB stars in a single starburst of 10^ single stars 
through one-dimensional linear interpolation of 
stellar masses. 

Figure [6] presents the dust masses produced by 
CE evolution in a single starburst of 10^ binary 
systems or by 10^ single stars. Compared with 
AGB stars, due to the high initial temperature 
and slow temperature decrease, the quantities of 
dust produced by CE evolution are negligible in 
the simulation with 7 = 0.3 at the base of the 
FGB. However, CEs produce significant amounts 
of dust or even dominate in the simulations with 
7 = 0.4 and 7 = 0.3 at the top of the FGB. 
This demonstrates the potential importance of CE 
evolution to the overall dust production in the 
ISM. 

Figure [7] shows how the dust masses produced 



by CEs or single AGB stars depend on the initial 
stellar masses. There are two peaks in the distri- 
bution of the dust quantities produced by AGB 
stars. The left peak is a direct consequence of 
the initial mass function that favours lower-mass 
stars, while the right peak originates from stars 
with masses of 3 — 4M (ti which are particularl y 
efficient dust producers (I Ferrarotti fc Gaill l2006h . 



The dust produced by CEs mainly comes from the 
stars with masses of 1 — 3 M© due to the initial 
mass function. 

7. Conclusions 

We have investigated the dust formation in CE 
ejecta constructing simple models for the evolution 
of the CE ejecta and using the AGBDUST code 
to simulate dust formation. The dust masses 
produced by CE ejecta greatly depend on the 
input parameters. Compared to the dust masses 
produced by AGB stars, they may be significant 
and could even dominate under certain conditions. 
This demonstrates that the contribution of dust 
produced by CEs to the overall dust produc- 
tion in the ISM needs to be taken into account. 
The progenitors of SNe la usually undergo CE 
evolution. Due to the generally expected long 
delay times of SNe la, it is difficult to observe 
the dust grains formed in CE ejecta via light 
echoes, but it may be possible to find light echoes 
produced by dust grains in SNe la in which their 
progenitors are composed of WDs and massive 
companions. Due to the short delay times, SNe lb 
may be good candidates for observing dust grains 
formed in CE ejecta. However, due to the still 
large uncertainties in modeling the CE phase, our 
conclusions are very preliminary and further work 
on its importance will be required. 
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